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Two loop mass effects in the static position space QCD-potential 
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The perturbatively calculable short distance QCD potential is known to two loops including the effect of 
massive quarks. Recently, a simple approximate solution in momentum space was utilized to obtain the potential 
in coordinate space. The latter is important in several respects. A comparison with non-perturbative lattice 
results is feasible in the overlap regime using light MS masses. This might be even more promising employing the 
concept of the force between the heavy color singlet sources, which can be easily derived from the potential. In 
addition, the better than two percent accuracy bottom mass determination from T-mesons is sensitive to massive 
charm loops at the two loop order. We summarize recent results using exact one loop functions and explicit 
decoupling parametrizations. 



In analogy to QED, the heavy quark potential 
Jljg] is of central interest in QCD as a measure of 
the strong coupling. In the long distance regime 
~ O(lfm) it contains the confining linear tail. 
With the inclusion of dynamical quarks on the 
lattice, one expects to enter the so called "string 
breaking" regime at distances of C(2fm) yielding 
the famous Yukawa potential H. At very short 
distances (< O.lfm) one recovers the Coulomb- 
part of the potential, modified by loop correc- 
tions like in QED. The overlap region is naturally 
of considerable interest since non-perturbative ef- 
fects might also enter [0 and possibly as a "nor- 
malization" for lattice calculations. The impor- 
tant point to notice here is that fermions with 
light MS-masses are easier to implement on the 
lattice. For a direct comparison, due to its smaller 
value, the force || between the static sources is 
also of special interest. For a recent review of lat- 
tice results on the forces in heavy bound states see 
Ref. H . From a phenomenological point of view, 
the importance of calculating mass effects in the 
perturbative part of the heavy quark potential is 
two fold. On the one hand, since we are dealing 
with a physical system, the corresponding phys- 
ical coupling definition (V-scheme) has several 
welcome advantages compared to massless (cal- 
culational) schemes. av{Q 2 /m 2 ) is an observable 
i.e. the flavor thresholds are analytic, display au- 



tomatic decoupling of heavy quarks and are inde- 
pendent of the renormalization scale to the order 
we are working. In addition, the physical scale 
of the problem is determined by the transfered 
momentum. Detailed discussions of the flavor 
threshold treatment using analyitc schemes are 
contained in Refs. ffl-ffl- On the other hand, 
the two loop mass corrections to the heavy quark 
potential are important for the better than two 
percent accuracy determination of the bottom 
mass [(To). In this case it is important not to 
treat the charm loops as massless but to con- 
sider the full massive calculation. Using the phys- 
ical T-meson for this purpose, the effect of the 
mass shift Srrib depends on (0i s |Vp(r, m c )|0i s ), 

where <f> ls = A- C""^" 3 ) 1 exp (-^-m b r) de- 
notes the Is ground state wave function of the 
T-meson and Vf the massive fermionic correc- 
tions to the potential. Note that in momentum 
space one would have an additional integration 
since each wave function depends on a different 
three momentum making the overall calculation 
prohibitive. 

We begin, by recalling the definition of the 
potential through the manifestly gauge invari- 
ant vacuum expectation value of the Wilson loop 
W T = {Q\TvVe^p(ig§ T dx^T a ) |0) of spatial 
extension r and large temporal extension T with 
gluon exchanges between the temporal lines. The 



path-ordering V is necessary due to the non- 
commutativity of the SU(3) generators T a . In 
the perturbative analysis through two loops con- 
sidered here, all spatial components of the gauge 
fields A£(r, ±T/2) can at most depend on a power 
of logT and are thus negligible here. Further- 
more, Wt — ► exp (— iTEo), where the ground 
state energy Eq is identified with the potential V. 
Thus we arrive at the definition: 



V(r, to) 



lim — 

T^ociT 



log(0|TrP exp fig I dx^A^A |0) (1) 

Writing the source term of the heavy charges, sep- 
arated at the distance r = |r — r'|, as 



JZ(x)=igv lt T a [6(x-r)-5(x-T>)] 



(2) 



and neglecting contributions connecting the spa- 
tial components, the perturbative potential is 
given by 



V(r, rn) = — lim — — 



log(0|TrTexp U d^A^J^x) 1 |0) (3) 

In the above equation v^ = J Mj o due to the purely 
timelike nature of the sources. The potential 
in momentum space is the Fourier transform of 
V(r). It can be calculated directly in momentum 
space from the on-shell heavy quark anti-quark 
scattering amplitude (divided by i) at the physical 
momentum transfer q, projected onto the color 
singlet sector. The potential can be used to define 
the effective charge ay ( the so-called V-scheme) 
through: 



V(Q,m) = -4vrC F 
V(r,m) = -C F 



av(Q,m) 

Q 2 

av(r,m) 



(4) 
(5) 



where Q = q = — q and both expressions 
above are related through a Fourier-transform, 
i.e. V{r,m) = J -§^V{Q,m) exp(iQr). At 
lowest order we obtain therefore the well known 
Coulomb potential Vc in each representation. 



In Ref. [fy] the two loop corrections in mo- 
mentum space were given in approximate form 
based on the reconstructed solutions of the Gell- 
Mann Low function with massive quarks obtained 
in Ref. || from the exact numerical results in Ref. 
|12| . In order to render decoupling of heavy fla- 
vors with pole mass m explicit one has to use the 
decoupling relation Q : 



<*<"»>(/*) = ai"'- 1 ^ 1 + - . "~' i log^- 






(mi 



- i] u 



7_ 

24 



19 M - 

— log — 
24 to 



1 , 2 F 
— log — 
36 TO 



(6) 



It is then useful to write the mass corrections to 
the potential in such a way that the light quark 
with mass m is included in the evolution of the 
strong coupling. Thus the mass corrections van- 



ish in the limits to 



and Q 2 



We 



find from the results of Ref. |11| the following 
expression: 

F NNL (Q, to) = 1/ NNL (Q, 0) + 6V N ^(Q, to) (7) 

where the first term on the r.h.s. is given in Ref. 
||14| correcting an error in the original calculation 
of the pure glue part of Ref. |jj|JL6fl . Denoting 
the one loop subtracted mass correction function 
by 



_ 5 
= 3 

with 

f(x) = V^-l 
we find 



1 & 

777/ 
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dxfjx) 
i Q 2 + Am 2 x 2 



2x 2 



2x 4 



(8) 



(9) 



SV NNL (Q,m) = V c (Q 




T F P 
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8 fm 2 \ /31 20 



^TlV 2 



c\ 



' dx 



16 
2Q 2 



m 2 \ 76 



x Q 2 + Am 2 x 2 
161 26 



•di 



' dx 



2Q 2 



114 19 



Cs 



log^ 

vnA 



x Q 2 + Am 2 x 2 
(10) 



The fitting constants are given by c 2 = 0.47±0.01, 
d-2 = 1.12±0.02 with ci, d\ fixed by the two condi- 
tions ci+di = 1 andcilog(4c|)+(ii loglog(4d|) = 
Hi + ffC3 i n order to ensure exact decoupling 
when using relation g. Eq. |o| can be used to 
obtain the corresponding results in coordinate 
space. Here we find analogously (with m = e 7 m): 



37T 



dxe~ 2mrx f(x) 



log(mr) + -+ 
6 



dxf(x) 



27i 



2^,2 



A) log 



4m x 

~ 2 



Fii(2mrx) — Ei(—2mrx)\ 
~fc A + ™T F n}j + A) y + 3 (log(mr) + |)x 
f/3 log(Ar) + — C X - -T F myJdxf{x)e~ 2mrx 

Ei(2nia:r)— 



\J x 1 



1 



— + x/(x) log - 

x^ x + V£ 2 - 1 



Ei(-2mxr) + log(4a; 2 ) ) -\ t- (log(mr) + - 

J 12 V 6 

57 /161 13, , ,„ , 
+ T I 228 + 19^ +1 °s( mr ) +ciEi(l,2c 2 mr) 



+diEi(l, 2d 2 mr) 



(11) 



where Ei(l,x) = J^ 00 exp(— i)^. The relation 
Ei(«) + Ei(-a:) = P / °°exp[(l - t)x]^ is also 
useful, with P denoting the principal value pre- 
scription. While Eq. [n] vanishes for m — ► 0, it 
does not for r — ► 0. It is a rather interesting fact 
that through the Fourier transform of Eq. [I(J a 
non-analyticnlinear mass term is generated which 
is furthermore enhanced by a factor ir 2 . It is im- 
portant to note that the linear mass term cannot 
be obtained by first expanding Eq. O and then 
integrating over the dispersion relation variable 
x, since it originates from large values of x. This 
means it originates from momenta smaller than 

1 Non-analytic refers here to the f act that Eq. Ila should 
be read as being proportional to Vm 2 since m is the only 
mass term entering the momentum space result in Eq. llOl 



the quark mass, i.e. it is of infrared origin. We 
find the following limit: 



r <l/m ^Cf f^(n,) 



(«S^)) : 



AV NNL (r,m) -m- 

31 „ 20 „ . /, ir 



1 



«£W 



2tt 



—C A - —T F ni + A log ^ - 41og(2) 



\ 4 76 

+3 1 + — (31 - 301og(2)) + — (ciC2 + di(fe) 

+C(m 2 r) (12) 

In terms of the MS-mass parameter m(/i), one 
only needs to use the relation 



m = m(fi) 



l + C F - 



7T V 2 



log 



M 



m(/i) 



(13) 



in the one loop term to obtain the NNL correc- 
tion. The size of the mass corrections relative to 
the Born Coulomb potential is displayed in Fig. 
[ljfor the charm quark mass m — 1.5GeV and the 
"natural" renormalization scale choice u = -K?, 
since the 7 terms are generated by the Fourier 
transform. In Ref. [O] it is shown that this scale 



choice is almost identical to the BLM-scale 18 
and thus consistently absorbs large renormaliza- 
tion group logarithms. The two loop running cou- 
pling is given by 



«T 



i(»)- 



4tt 



/ 



A log 



^QCD \ 



A lp g ( lQ g A 
A) log 



(14) 



A 2 



where we normalize the QCD-scale parameter 
h-QCD such that a-^(Mz) = 0.12 which corre- 
sponds to Aqcd = 0.25GeV and we keep m = 4 
fixed. The first two terms of the /3-function are 
gauge invariant and scheme independent in mass- 
less renormalization schemes and are given by 
A = fC A -±T F rn and A = f C\ -f C A T F m- 
4C F T F ni. In QCD we have C A = 3, C F = | and 
T F — i . The effect is at the several percent level 
and increases above 10 % for distances of 0.1 fm 
(~ 0.5 GeV -1 ). It vanishes for small r since we 
are displaying the corrections relative to Vc- In 
absolute terms we checked numerically that Eq. 
ft2| is reproduced for small distances by the full 




Figure 1. The size of the charm-mass corrections 
in Eq. pll relative to the Coulomb potential. The 
choice of the natural renormalization scale is in- 
dicated in the figure. The distance r between the 
sources is in GeV -1 . 



result in Eq. [Tl| The inclusion of these charm 
quark corrections into a full sum rule determina- 
tion of the bottom quark mass in Ref. [[19[ led to a 
shift of -30 MeV and yields rn b (m b ) = 4.17±0.05 
GeV, which in light of the small error is a signif- 
icant contribution. 

Instead of considering the potential, or 
av(r,m), it is also of interest to study the cou- 



„2 d(ay(r,m)/r) 
dr 



which is de- 



pling ap(r, m) = 
fined from the force between the static quarks. 
In general aF(r^m) is smaller than ay(r, m) or 
even «5Ss(/i = 1/r) |ll[, which makes it suit- 
able as an expansion parameter in perturbative 
calculations. It is also useful in lattice calcula- 
tions to determine a s |2(j,|2l]]. The force is given 
by F(r,m) = — Lpml anc [ the massless case is 
taken from Ref. [fill : 



F NNL (r,0) = -Cf 



ff&oL , «&>(/*) 



>o 



2/3 log(Ar)+/i + 



?M 



■in 



All 



4/3 2 log 2 (Ar) 



ft + 2A)/i 



log(pr) + f 2 



(15) 



where 

fl = ~~9 A + 9 FUl 



(16) 
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16C 3 C F T F m 
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"8T~ 27 5 



7T 2 (T F n,) (17) 



From Eq. O we thus find for the mass corrections 
the following expression: 

<5F NNL (r,m) = - dl ° g ^ c ^ SV^^m) - V c \ 
f(x)\2mx(p 



~y Ca + t pni 





-ih 



log(mr) + 



3 / 31 10 

-- f A)log(Ar) + — CU - y T ^ 



dxe 



-2mra 



/(x) 2ma; log(4x 2 ) - Ei(2mrx) 



Ei(— 2mrx) H — - + a;/(x) log 



+ \Jx 2 - 1 



+ e' 



c ) j + 1 (log(mr) + 



57 fi 



2comr j — 2/i27nr"\ 

— die l 



(18) 



The size of the charm-mass corrections relative to 
the Born term is presented in Fig. [| for the same 
scale choice as in Fig. [j] It can be seen that 



(i=1/r/e T 
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Figure 2. The size of the charm-mass corrections 
in Eq. O relative to the lowest order force from 
the Coulomb potential. The choice of the natural 
renormalization scale is indicated in the figure. 
The distance r between the sources is in GeV -1 . 



the effect is less than half of that for the mass 
corrections to the potential and start to increase 
more rapidly at distances over 0.1 fm. 

In summary, we have calculated the massive 
quark corrections to the static QCD potential in 
coordinate space at the two loop level. The re- 
sults presented here use the exact one loop vac- 
uum polarization functions and a dispersion rela- 
tion fit based on the results of Ref . JO] . The un- 
certainty is estimated at the percentile level from 
comparisons with the exact Monte Carlo results 
in momentum space of Ref. p2]. For the bot- 
tom mass determination, the inclusion of a mas- 
sive charm quark in the T potential is significant 
due to two reasons. One factor is that the effec- 
tive physical scale depends parametrically on the 
charm mass, leading to a large value of the cou- 
pling. The other reason originates from the fact 
that in the relation between the potential contri- 
bution to the static energy and the pole mass, 
there is an uncanceled non-analytic linear mass 
term [ |10|JlS| ], whose origin is the Fourier trans- 
form. Together, these two effects lead to a shift 



of -30 MeV and M b (m b ) = 4.17 ± 0.05 GeV. In 
addition we have calculated the mass corrections 
for the force between two static color charges in 
a singlet state at the two loop level. The size 
of the effect in general is smaller than for the 
potential but still significant at larger distances. 
At small distances, the linear r-independent mass 
term from the potential drops out and thus leads 
to smaller corrections in this regime. 
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